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Abstract 



The ALEPH data on the (non-strange) vector and axial- vector spectral functions, ex- 
tracted from tau-lepton decays, is used in order to search for evidence for a dimension- 
two contribution, C2 v,A, to the Operator Product Expansion (other than d = 2 
quark mass terms). This is done by means of a dimension- two Finite Energy Sum 
Rule, which relates QCD to the experimental hadronic information. The average 
C2 = {Civ + C2a)/2 is remarkably stable against variations in the continuum thresh- 
old, but depends rather strongly on Aqcd- Given the current wide spread in the 
values of Aqcd, as extracted from different experiments, we conservatively conclude 
from our analysis that C2 is consistent with zero. 



*Work supported in part by the Volkswagen Foundation 



The Operator Product Expansion (OPE), extended beyond perturbation theory, is one 
of the pillars of the successful QCD sum rule method used extensively to link QCD and 
hadronic physics |T|. In analyzing e.g. two-point functions, one calculates QCD perturba- 
tive contributions up to a desired order in the running strong coupling, and then includes 
non-perturbative effects parametrized in terms of a series of vacuum to vacuum matrix 
elements of local, gauge invariant operators built from the quark and gluon fields entering 
the QCD Lagrangian. These so called vacuum condensates encompass the long distance 
dynamics. They are multiplied by Wilson coefficients, calculable in perturbation theory, 
which contain the short distance information. For a given dimension, these terms fall off as 
powers of the (Euclidean) momentum transfer q 2 (q 2 < 0). The lowest naive dimension in 
QCD is d = 4, corresponding to the gluon condensate and to the product of quark masses 
and the quark condensate. In the standard OPE approach there are no dimension-two cor- 
rections, other than the well known quark mass insertion contributions of the form m 2 /q 2 . 
Except possibly in the strange-quark sector, these terms can be safely neglected for light 
quark current correlators. On the other hand, some specific dynamical mechanisms have 
been suggested as potential sources of d = 2 corrections to the OPE, e.g. renormalons 
0] or a tachyonic gluon mass Q. If present, these dimension-two corrections are not ex- 
pected to have much impact on most of the existing phenomenological predictions from 
QCD sum rules, as the level of precision is seldom better than 10-20 %. However, they 
may have a non-negligible impact on the extraction of the QCD strong coupling at the 
tau-lepton mass scale, a s (M 2 ). In fact, current claims from this precision determination 
|| have been questioned []5| on the grounds that they rely on the assumption of no d = 2 
corrections to the OPE (other than quark mass insertions). Attempts have been made to 
determine phenomenologically the size of potential dimension-two terms using information 
on (i) the vector spectral function as obtained from e + e~ data ||, and (ii) the vector 
and the axial- vector spectral functions extracted from tau-lepton hadronic decays [0, as 
measured by the ARGUS collaboration ||. The analysis based on tau decay data has the 
advantage of relying on two independent quantities (d — 2 corrections are expected to be 
chiral-symmetric), which can be constrained further by using theoretical information such 
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as the first Weinberg sum rule. The conclusion from this determination [7|] was that the 
ARGUS data supported a non-zero dimension-two term in the OPE, which was consistent 
with a dependence on the QCD scale of the form: C 2 oc Aq CD . With the advent of the 
ALEPH precision data || on semileptonic tau decays it is now possible to analyze a vari- 
ety of theoretical issues involving the vector and axial-vector spectral functions. One such 
issue is, precisely, the possible presence of a dimension-two contribution in the OPE, which 
is the subject of this note. These spectral functions are related to the discontinuities in 
the complex energy plane of the two-point functions involving the vector and axial-vector 
currents 

n^(g 2 ) = ^ / d 4 xe^ <o|T(v» K+(o))io > 

= {-9p,u q 2 + q»qv) ny(g 2 ) , (1) 

n£V)=i / d'xe^ <0\T(A,(x) 4(0))|0> 

= (-9^ q 2 + q^qu) iWg 2 ) - q^q v n (g 2 ) , (2) 

where =: (uj^u — dj^d) : /2, and =: {u^/y^u — d^^^d) : /2. Considering these 
(charge neutral) currents implies the normalization Im ITy = Im Ua = 1/87T, at leading 
order in perturbative QCD. In order to determine the size of a potential dimension-two 
contribution to the OPE we consider the following Finite Energy Sum Rule (FESR) |10| 



t - 871 f t \ a ° 2V > A , TP ( \ (o\ 
hv,A = / PvA s )ds = \- F 2 {so) , (3) 

s Jo s 
where Ci v,a is the potential d — 2 term, so is the continuum threshold, and F 2 (so) is the 
radiative correction, identical in the vector and axial-vector channels, which is obtained 
after a straightforward integration of the perturbative QCD results of [TT|; to four- loop 
order it is given by 

« W -i + ^ + (^)>-^-f) 

+ i =^ } | (ln2i _ lni _ 1) + |_ 2 ( F 3-f)01n i + F 4 -^-f + f], (4) 



with 

where L = ln(s /A^ CD ), and for three flavours: ft = -9/2, ft = -8, ft = -3863/192, 
F 3 = 1.6398, F 4 = —10.2839. In writing Eq. (4) use has been made of the result [|ll|] 

«i 3) (*o) = ggj(fg) + (o^fo) \ 2 / z^2 lnL \ 

7T 7T I 7T j I ft i 



+ (^Jli (ln2L - lnL - 1)+ fJ' (6) 

and an expansion in powers of a^' is to be understood. Alternatively, one may choose not 
to expand in such a way; in this case the radiative correction becomes 

F2(3o) . 1 + ^M + (^M) 2 ^_|) 

where a^ 3 ^(so) is given by Eq. (6). Numerically, these two alternatives have a non- negligible 
impact on the final result for C2, as will be discussed later. 

The quark mass insertion term, which contributes to this dimension-two FESR is of the 
form 



,„,2 .;,2 

C*2m = — 3— j- 



mi + m l dj 



Using current values of the up- and down-quark masses, this term is negligible. 



It should be stressed at this point that the FESR are ideally suited, in principle, to extract 
the values of power corrections of a given dimensionality. Ignoring gluonic corrections to 
the condensates, the FESR involving py.A with kernel s N (N=0,l,2,...) project out only 
condensates of dimension d = 2,4,6, ... .In other words, in the FESR of lowest dimension 
all condensates of d = 4,6,... decouple. This should be contrasted with e.g. Laplace or 
Gaussian sum rules which receive contributions from all possible condensates. Since the 
numerical values of these power corrections are not well known, these other sum rules in- 
troduce an unnecessary additional uncertainty. On the other hand, the fact that FESR 
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tend to emphasize the high energy region, where the ALEPH data have larger errors, is of 
no importance here, as the uncertainty in C 2 turns out to be overwhelmingly dominated by 
the uncertainty in Aq C d, and to a lesser extent, by the way the perturbative expansion is 
organized. The experimental error in the hadronic integral in Eq.(3) can be safely neglected. 



We show now that the dimension- two terms obtained from Eq.(3), i.e. 

2 f s ° 

C 2 v,A = 8n / Pv,A( s ) ds ~ s o-F20o) , (9) 
Jo 

are actually identical in the vector and axial-vector channels, provided one takes the chiral 
limit. This result would follow trivially from e.g. the first Weinberg sum rule, provided 
this sum rule would be saturated by the data for s < 00 (actually, s < Mj? in the case 
of tau decay data). However, this is not the discussed in flT2"| . Instead, the data 



saturate much better the modified sum rule 12 



r s ° s 

/ (1--) [ Pv (s )- PA {s )]ds = 0, (10) 
Jo s 

Here, pa already contains the pion pole, i.e. 

p A {s) = fj{s) + p A {s)\ RES (11) 

where Pa(s)\res is the resonance part of the spectral function. We make use of Eq.(9) in 
Eq.(10), and invoke the dimension-four FESR 



hv,A = — I pvA s ) sds = — o 2 ( 12 ) 

Sq Jo £ Sq 



where both the radiative correction F 4 (s ), and the dimension- four condensate (equal to the 
gluon condensate in the chiral limit), are identical in the vector and axial- vector channel. 
One then obtains 



f ° (1 - -) {p v (s) - PA (s)]ds = = -^(CW - C 2A ) , (13) 

JO Sq 07T Z 

which completes the proof. 
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We proceed now to use the ALEPH data in Eq.(9) in order to check for any evidence for 
a dimension-two operator. The results of a numerical evaluation of the r.h.s. of Eq.(9), 
using the expanded form of the QCD integral, Eq. (4), and a fit to the ALEPH data as 
described in |L2|], are plotted as a function of s in Fig.l, for Aq C d = 300 MeV. Curve 
(a) is the result in the vector channel, curve (b) in the axial-vector channel, and curve 
(c) is the average between the two, i.e. C 2 = (C 2 y + C 2 ^)/2. It can be seen that the 
results for Civ and C 2 a, considered individually, are rather unstable against variations in 
So- This behaviour is simple to understand. When integrating the data, the vector integral 
at So — 1 GeV 2 has already picked up the contribution from the (narrow-width) rho-meson, 
while the axial- vector spectral function is still relatively small there. The contribution from 
the (broad-width) ax-meson is important only for so > 1.5 — 1.6 GeV 2 . For this reason, 
the hadronic integral approaches the theoretical or QCD integral from above in the vector 
channel, and from below in the axial- vector channel at so < 1.5 GeV 2 . The expectation 
C2V = C2A is only true asymptotically. However, making use of this constraint, it is 
natural to consider instead the average value. This turns out to be remarkably stable, 
and allows for a precise determination of Ci- Perhaps this should not be surprising, as 
the dimension-six four-quark condensate contributes with different signs to the vector and 
axial- vector correlators, and there is a tendency to an overall cancellation between the sum 
of the gluon and the four-quark condensates in + (which, however, enter the sum 
rule under consideration only via radiative corrections). Nevertheless, there is a strong 
dependence of C 2 on Aq C d as discussed next. The size of the current error bars in Aq C d 
makes it the dominant source of uncertainty in the determination of C 2 ; in comparison, the 
small experimental uncertainty in the hadronic spectral functions plays a negligible role. 

First, we have found a strong dependence of C 2 on the order at which the radiative cor- 
rection F2, Eq.(4), is computed. As the right hand side of Eq.(9) is the (small) difference 
of two similar numbers, it is not surprising that the transition from one-loop to two-loops 
of perturbative QCD in Eq.(4) leads to relatively large differences. In other words, the 
result for C 2 obtained using Eq.(4) truncated to one loop, i.e. using only the first term in 
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that equation, differs substantially from the result for C2 if the truncation is to two-loops, 
i.e. including the first two terms in Eq.(4). Higher order corrections (three-, four-, and 
five-loops) show, however, the expected convergence, i.e. there is little difference between 
the values of C2 obtained after truncation at the three-loop, four-loop, and five-loop level. 
In the latter case, Eq.(7) becomes 

2 / 



^2 (so) 



1 + qW( ao ) + M 4 W 



+ 



a 



(4) 1 



7T 



SO J 



7T 



7T 



+ 



a 



(4) 1 



so) 



1 r 



7T 



PI 



71 
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(14) 



where a^(s ) is given by [|13 



(3), 



7T 



and 



h = — 



44 



7T 



149753 



^(lnL--lnL-21nL + -) + 3- r lnL + - 



(15) 
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f3564C 3 



1G7S3G1 G.-.0S .50065 6472^ ^ 2 t 1093 m3 



162 



27 



162 



-C 3 )<+ 



729 



(16) 



and £3 = 1.202. While the QCD beta-function is known to four loops |T3[], the imaginary 
part of the vector (axial-vector) correlator to five-loop order involves the unknown constant 
k 3 in Eq. (14). We have estimated this constant assuming a geometric series behaviour for 
those constants not determined by the renormalization group, i.e. k 3 ~ k\jk\ ~ 25, with 
k\ = F 3 , and /C2 — F4 + 7T 2 (3\j '12. This is in good qualitative agreement with other estimates 



i~4| . In this case we find that the difference between the four-loop and the five- loop results 



for C2 is again quite small. Nonetheless, the four- and five-loop determinations (represented 
in Fig. 2) depend strongly on the QCD scale parameter Aqcd (here Aqcd = A.jj^(n F = 3)). 
As Aqcd is varied in the range JT^| 300 MeV < Aq C d < 400 MeV, the extracted value 
of C2 remains small and changes sign after crossing zero for Aqcd — 330 — 360 MeV, for 
typical values of sq in the stability region (so — 1.5 - 2.7 GeV 2 ). In Fig. 3 we show the 
dependence of C 2 on Aq C d for a typical value s = 2 GeV 2 at the center of the stability 



region. Such a strong dependence of C 2 on Aqcd could be a welcome feature, if for some 
theoretical reason one could rule out completly the presence of a dimension-two opera- 
tor in QCD. One could then determine Aqcd with great accuracy. Finally, the difference 
in the QCD integral after expanding in powers of ajp, as in Eq.(4), or not expanding 
them, as in Eq.(7), though small, does affect the value of C 2 , as this is the result of the 
subtraction of two very similar numbers. For Aqcd as determined by the PDG |I5fl , i.e. 
Aqcd = 389 ± 35 MeV, we find C 2 = -(0.08 ± 0.06) GeV 2 if the expansion in a« is done 
as m Eq. (4), and C 2 = -0.05 ± 0.05 GeV 2 if one does not expand, as in Eq.(7). Alterna- 
tively, we consider the value of a s extracted in JJ, which implies Aq C d = 367 ± 40 MeV, 
for three flavours. This analysis was performed assuming C 2 = 0. In view of our re- 
sults, this assumption is not justified a priori. However, such a value of Aqcd leads to 
C 2 = -(0.05 ± 0.06) GeV 2 using Eq.(4), and C 2 = -0.03 ± 0.06 GeV 2 using Eq.(7), which 
are consistent with zero, and thus justifies the assumption a posteriori. All the above nu- 
merical calculations were performed using an analytical chi-squared fit to the ALEPH data, 
and also by direct integration of these data , as they are available in bin-interval format 
||. The difference between the hadronic integrals in Eq.(3) from these two procedures is 
less than 1 %, which provides a reasonable estimate of the error in integrating the data. 
While this small difference translates into a larger difference in the values of C 2 , due to 
the cancellation between the hadronic and the QCD integrals, the conclusion remains the 
same, i.e. that given the present uncertainty in Aqcd the ALEPH data imply a value of 
C 2 consistent with zero. This conclusion is rather different from that reached in a previous 
analysis J7J, along the same lines as here but using instead the ARGUS tau-decay data || 
to determine the hadronic spectral functions. This discrepancy is due in part to the much 
larger error bars of the ARGUS data. However, this does not fully explain the differences, 
which are mostly due to different values of the vector and axial-vector spectral functions 
at, and in the vicinity of the respective resonances (rho- and ai-mesons). This translates 
into different areas under the hadronic spectral functions. For instance, the ARGUS data 
saturated the first and second Weinberg sum rules reasonably well |ToT|, while this is no 



longer the case for the ALEPH data||, [12 
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Figure Captions 

Figure 1. The right hand side of Eq.(9) as a function of the continuum threshold s 
in the vector channel (curve (a)), the axial-vector channel (curve(b)), and the average 
C 2 = (C 2 y + C 2 a)/2 (curve (c)), all for Aq C d = 300 MeV. The expanded expression, Eq. 
(4), has been used. 

Figure 2. The average C 2 as a function of the continuum threshold so for Aqcd = 300 MeV 
(curve (a)), and Aq C d = 400 MeV (curve (b)). The expanded expression, Eq.(4), has been 
used. 

Figure 3. The dependence of the average C 2 on Aq C d for a typical value s = 2 GeV 2 at 
the center of the stability region. The expanded expression, Eq.(4) has been used. 
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Figure 1: 
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Figure 2: 
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Figure 3: 



13 



